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TECHNICAL NOTE NO. 17^7 


DETERMINATION OF COUPLED AND UNCOUPLED MODES AND 
FREQUENCIES OF NATURAL VIBRATION OF SWEPT AND 
UNSWEPT WINGS FROM UNIFORM 
CANTILEVER MODES 

By Roger A. Anderson and John C. Houbolt 


SUMMARY 


A solution is presented for the coupled and uncoupled inodes and 
frequencies of wings (swept or unkwept) mounted on a fuselage. The 
deflection and twist of the wing is expanded in terms of the modes of a 
uniform cantilever team and the energy method is used to derive the 
characteristic equations describing symmetrical and antisymmetrical 
modes of vitiations. Application of these equations to various types of 
wing vibration is illustrated ty numerical examples . The numerical 
examples chosen were susceptible to exact solution, and the results show 
that accurate modes and frequencies can te obtained ty a method in which 
low -order determinants are used. 


INTRODUCTION 


Except for certain Idealized cases, the natural vibration modes and 
frequencies of airplane wings (swept or unswept) cannot te found by 
exact analysis, and thus approximate methods of solution must te used. 
Such a solution is presented for the general problem of coupled tending 
and torsional vibration of a nonunlform wing mounted at an angle of 
sweep on a fuselage. The energy method is used to derive two sets of 
linear characteristic equations, one for sy mme trical modes and the other 
for antisymmetrical modes. These same equations also lead to solutions 
for the coupled and uncoupled modes and frequencies of the unswept wing 
and fuselage. The equations are solved ty the Crout method which per- 
mits Independent calculation of any desired mode and frequency. 

In this solution the deflection and twist of the wings are assumed 
to conform to elementary team theory. Such an assumption may or may not 
lead to an analysis which is applicable to wings having low aspect 
ratio, especially when appreciable sweep is present, because the effects 
of the distortions in the vicinity of the root of the wings are not 
fully understood. The present analysis considers the fuselage to te a 
rigid tody, tut the analysis may te extended to treat an elastic 
fuselage. 
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Tlie important feature of the method presented herein is the simpli- 
fication that results from use of the natural modes of a uniform canti- 
lever beam as expansions for the deflection and twist of the vibrating 
nonuniform free -free wing. With the addition of appropriate rigid -body- 
displacement terms, these expressions may be made to satisfy all the 
geometrical boundary conditions for both symmetrical and antisymmetrical 
wing vibration. Analyzing idealized structures for which exact solutions 
could be made shows that only a few terms are needed in the expansions 
to obtain good accuracy. 


SYMBOIS 


L length of semi span along elastic axis 

E Young 's modulus of elasticity 

G modulus of elasticity In shear 

I bending moment of inertia of wing cross section perpendicular 

to elastic axis 

Ija polar moment of inertia of mass per unit length of wing about 

elastic axis 

J torsion constant for cross section perpendicular to elastic 

axis of wing 

Ipp one -half of pitching polar moment of inertia of fuselage mass 

about elastic axis of wing at root 

Ipp one -half of rolling polar moment of inertia of fuselage mass 

n about fuselage longitudinal axis 

x distance along elastic axis of wing, measured from center 

line of fuselage or root of wing 

y deflection of elastic axis of wing with respect to its static 

equilibrium position} positive upward 

<p angle of twist of wing cross section with reference to. its 

static equilibrium position} positive in st allin g direction 

9 angle of pitch of fuselage 

$ angle of roll of fuselage 

angle of sweep, measured between wing elastic axis and line 
perpendicular to fuselage longitudinal axis 


A 
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m mass of wing per unit length 

TOy total mass of wing 


mg* 

*m 

co 


one -half of mass of fuselage 


ratio of fuselage mass to vi mg mass 




) 


circular frequency of natural mode of vibration, radians per 
second 


p n circular frequency of nth natural tending mode of uniform 

cantilever team. 

q n circular frequency of nth natural torsional mode of uniform 

cantilever team 

f frequency of natural mode of vibration, cycleB per second 

e distance between mass center of wing crosB section or points 

of mass concentration and elastic axis of wing} positive 
when mass center lies forward of elastic axis 

ep distance between mass center, of fuselage and elastic axis of 

wing at root} positive when mass center lies forward of 
elastic axis 

e a distance between point of spring suspension and elastic axis 

of wing at root} positive when point of suspension lies 
forward of elastic axis 


a elastic spring constant 

a n coefficient of nth term In expansion for y 

bja coefficient of mth term in expansion for 9 


i, m, n 

TnW 

y n "W 

<P n ' (*) 


integers (1, 2, 3, . . .) 

nth bending mode of uniform cantilever beam 
second derivative with respect to x of y n (x) 
nth torsional mode of uniform cantilever beam 
first derivative with respect to x of <p n (x) 
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ENERGY expressions and nefeection functions 


In an analysis to determine the inodes of Titration of wings it is 
usually sufficient to consider the equilibrium of one -half the fuselage 
apd the wing semispan only. The airplane is divided along its longi- 
tudinal axis with a coord in ate system assigned as shown in the following 
configuration : 



In this analysis the fuselage' is assumed to he rigid. Its motion is 
therefore represented hy rigid-hody displacement and rotation. 

For Titration of this system, the energies considered are the 
tending, twisting, and kinetic energies of the wing semispan and one- 
h»i-f the kinetic energy of the fuselage. At maximum displacement of the 
wing the sum of the strain energy of tending and twisting is given hy 
the well-known expression 



dx 


( 1 ) 
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The kinetic energy of the -wing as it passes through the equilibrium 
position is given by the equation (see appendix A for derivation) 



1 L 

PL 

< 

H 

II 

my 2 dx + as 2 

,.,2 

me<py dx + — 

J 

0 . i 

o 2 t 


PL 


VP 2 dx 


( 2 ) 


0 


One-half the kinetic energy of the fuselage is given similarly hy the 
expression (see appendix A for derivation) 



(3) 


It can he shown hy geometry that the pitching angle 0 and the rolling 
angle $ are related to the angle of twist of the wing at the root cross 
section, the slope of the wing at the root, and the angle of sweep hy the 
following relations 


0 


^cp cos A 


dy 

dx 


sin A 


'x=0 


(4) 


$ 


q) sin A + 


dy 

cos A 

dx 


x=0 


(5) 


The conditions leading to these geometrical relations were pointed out 
hy Lawrence in reference 1- 

By the energy method, expansions are chosen to represent the 
deflection and twist of the wing. In accordance with a principle of 
this method, which apparently is not very well known, the expansions 
need only satisfy the geometrical houndary condltipms of the problem. 

The dynamical boundary conditions (shear, moment, torque) are auto- 
matically satisfied in the minimization process. A brief account of 
this feature is given in reference 2. 

Convenient expansions for the deflection and twist are the following 
equations : 
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y(z) = Eq + a^, J + aqyq + + • • • + &n7 n + 


q)(z) = "bo + bqcpq + b^qjg + * • • ’V'Pni + 


( 6 ) 

(7) 


in which the coefficients a n and bm may he taken as generalized 
coordinates. The coefficients a^,, a^, and h Q are included to allow 
the wing to have a rigid -body deflection, a linearly varying deflection, 
and a rigid-body rotation. The y n 's and ij^'s are bending and 

torsional modes, respectively, of a uniform, cantilever beam written in 
terms of a unit tip deflection or unit tip rotation. Some of the reasons 
for choosing these modes will be discussed in the section entitled 
'discussion of Method." With the use of the expansions ( 6 ) and ( 7 ) the 
geometrical boundary conditions at the wing root are determined solely 
by the values of the coefficients a Q , a Q , and b Q . (.The coefficient a Q 

determines the wing root deflection, a Q determines the slope of the wing 

at the root, and. b Q determines the angle of rotation of the wing root.) 
The simplification resulting from this method of dealing with root 
boundary conditions will be made evident in the section to follow. 

Upon substitution of equations (4), ( 5 ), ( 6 ), and ( 7 ) into 
equations (1), ( 2 ), and ( 3 ), the energies U, Vq, and Tg are expressed 

in terms of the unknown coefficients b^ and the unknown frequency <n. 

(The coefficients a Q , a Q , and b Q do n<5t appear in the expression 
for U.) 

For the case of a flexible fuselage, expansions similar to 
equations ( 6 ) and ( 7 ) can be made for the deflection or twist of the 
forward and rearward sections of the fuselage. These expansions are then 
used with appropriate energy terms similar to equations (l) and ( 2 ). 


BOUNDARY CONDITIONS AND CHARACTERISTIC 
EQUATIONS OF SWEPT WINGS 


Because of the structural symmetry existing about the longitudinal 
axis of the fuselage, the structure vibrates in natural modes which can 
be classed as either symmetrical or anti symmetrical about this axis. 

For convenience, the boundary conditions defining each type of vibration 
will now be introduced. 

For symmetrical vibration, the boundary condition at the wing root 
is that the fuselage shall not roll, or 
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4 = 


<p sin A + ^2. cos 
dx 


= 0 


/x=0 


(8) 


Use of equations (’ 6 ) and ( 7 ) permits this boundary condition to "be 
expressed "by the following simple relation "between the coefficients a Q 
and "b 0 


a,-, = -t 0 L tan A 


(9) 


Elimination of Sq from Yq and Yg "by means of this relation leads 
to the solution for symmetrical modes and frequencies. 

O 

For antisymmetries!, vibration, the constra ini ng relations at the 
wing root are that the deflection is zero and that the fuselage shall 
not pitch, or 


(y) 


x=0 


= 0 


( 10 ) 


and 


0 = | qp cos A - — sin A] =0 

V ^ Jx=0 


( 11 ) 


By use of equations (6) and (7 ) , the "boundary conditions may "be written 


a 0 = 0 


( 12 ) 


^o 


— tan A 
L 


(13) 


Substitution of these relations into the expressions for Yq and Y 2 
leads to the solution for antisymmetrical modes and frequencies. 

The characteristic equations for modes and frequencies of Titration 
(symmetrical and antisymmetrical) can he found ty minimization of the 
expression Yq + Yg - U with respect to the unknown coefficients aq 

and tq. For example, the equation obtained ty minim ization of the 
energy expression with respect to Sq is 
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^( 1 % + m^) + a^Aoi + &2AQ2 + • • • a n Aon 


+ b Q 



+ 


jaqt.m 

cos A 1 


tan A 


+ 


"biFoi + 1 02 ? 02 . + • • * t>mF0m = 0 


(lb) 


The sets of linear homogeneous equations derived in this way are shown 
in a somewhat condensed form, for symmetrical and antisymmetrical vibra- 
tion in tables I and II. The constants Ain. B in. C in. • • • shown in 
these equations represent integrals which are defined in appendix B. 

Nor any problem, the limi ting values of the subscripts n and m of 
a^ and are chosen in accordance with the number of modes being 

represented J>y the expansions. The examples to be presented subsequently 
herein will serve as a guide in this choice. The next section will show 
that the equations which characterize coupled and uncoupled vibration of 
an airplane with unswept wings may be obtained s imply by crossing out 
certain rowB and columns in table I or table II. 


CHARACTERISTIC EQUATIONS NOR VIBRATION 
OF UNSWEPT "WINGS 


The equations for vibrating unswept wings may be obtained directly 
from tables I and H simply by letting the angle of sweep A equal 
zero. The equations for the various types of vibration - coupled 
bending and torsion, uncoupled bending, and uncoupled torsion - may 
then be found as described in the following sections. In all these 
cases it is understood that A = 0. 


Coupled Bending and Torsion 

Symmetrical and antisymmetrical free -free vibration . - The equations 
obtained for A = 0 describe symmetrical and anti symmetrical free -free 
vibration of an unswept wing and fuselage. 

Vibration as a cantilever . - The deflection, slope, and twist at the 
wing root are zero and therefore ao = a^, = b Q = 0. The characteristic 
equations for this type of vibration may be found from table II by 
crossing out the column headed by a^> and the row for Sq. 


Uncoupled Bending Vibration 

If wing twist and fuselage pitching are prevented, the twisting 
component <p is zero and consequently the "fc^'s must be zero. In 
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tables I and II, therefore, the upper left-hand quadrant, -which involves 
only the a^s, gives the equations for uncoupled bending vibration. 

The several types of uncoupled bending vibration may be found by dropping 
out certain of the remaining equations as follows: 

Symmetrical free-free bending .- For symmetrical free -free bending 
the slope at the root is zero and therefore aj-, = 0. The upper-left 
quadrant of table I yields the desired equations. 

Antisymmetrlcal free -free bending .- For antisymmetrical free -free 
bending the deflection at the root is zero, and therefore elq = 0. When 

A = 0 t he upper-left quadrant of table H yields the desired equations. 

Cantilever bending . - For cantilever bending the' deflection and 
. slope at the root are zero and therefore Sq = a^ = 0 . The equations 
for this type of vibration may be found from the upper-left quadrant of 
table I by crossing out the column headed by a Q and the row for a Q . 


Uncoupled Torsional Yibration 

If wing deflection* and fuselage translation are prevented, the 
deflection component y is zero and consequently the a n 's must be 
zero. In tables I and II, therefore, all the rowB (equations) derived 
from the minimization of the ap 's and the column headed by a n 's may 
be crossed out. The remaining equations, the lower-right quadrant, 
then define uncoupled torsional vibration of a wing-fuselage system. 

Symmetrical free -free torsion . - For symmetrical free -free torsion 
all the equations in the lower -right quadrant of table I are used 
(A=0). 

Cantilever or antisymmetrlcal free -free .torsion . - For cantilever or 
antisymmetrical free -free torsion the twist at the root is zero and 
therefore b Q = 0. The equations for this case are found in the lower- 
right quadrant of table II . 

SOLUTION OF CHARACTERISTIC EQUATIONS 


It has been shown that the equations characterizing various types 
of wing vibration may be found in a simple manner from the general sets 
of equations shown in tables I and H. A solution to a particular set 
of equations consists in determining values of cd, a^ and which 

define the natural frequencies and modes. 

Since the equations are homogeneous , values of a^ and b^ other 
t han zero which satisfy the set of equations can be found only when the 
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determinant of the set of equations is zero. The determ i nant, -which is 
simply the "bracketed terms in tables I and H, contains the unknown 
frequency coj the values of co which cause the deter minant to "be zero 
are the natural frequencies of vibration. In order to .deter mine the mode 
associated with a given natural frequency, one of the unknown coeffi- 
cients a n or "bn is given an arbitrary value (unity, for convenience) . 
Mathematical consideration indicates that one of the equations can now 
"be discarded, "but it- has "been found "best simply to add two of the equa- 
tions together and then to solve the resulting set of nonhomogeneous 
equations simultaneously to obtain the relative values of the other coef- 
ficients. The value of frequency used in the solution of these equations 
is, of course, the frequency value which caused the determinant to be 
zero. With the coefficients known, the mode associated with the given 
value of frequency is obtained directly from equations ( 6 ) and ( 7 ) . 

The values of o> satisfying the frequency determinant may be found 
by several methods. Perhaps the simplest way to locate a frequency root 
is to evaluate the determinant for a number of trial values of 00 in 
the expected vicinity of a natural frequency and to plot a curve of o> 
against the value of the determinant. In most cases, the value of cn 
giving a zero determinant can be obtained from the results of three or 
four evaluations. The evaluations may be performed by the Crout method 
of solving determinants. (See reference 3*) The Crout method yields 
solutions rapidly and provides for a running check which minimizes the 
possibility of computational error. With the procedure Just outlined, 
any desired frequency root and mode can be found independently of the 
other frequencies and modes. 


COMPUTATION OF CONSTANTS 


Before the characteristic equations for vibration of the wing can 
be solved, the constants A 1riJ> B^, C^, . . . appearing in these 

equations and defined in appendix: B must be determined from the physical 
properties of the wing. Appropriate numerical and graphical integration 
procedures are known for the calculation of such integrals. In order to 
use the procedures, it is necessary to know the physical properties m, 

me, Ijh, El, and GJ at a numb er of representative stations 2L along 

Xj 

the wing and also the numerical values of the modal functions y n (x) 
and (^(x) and their derivatives y^" and at these stations. 

For convenience, y n (x), y n "(x), %(x), and q^'(x) (n = 1, 2 , ... 5) 


have been computed at ten stations 



0.1 




0*2 


7 


0.9, 1.0 


) 


and 


are presented in tables ITT and IV. The values at stations other than 
those listed in the tables may be determined from plotting the modal 
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functions and reading off the desired values. The constants A-j-n, B^, 
C iri , . . . are then found "by multiplying the physical constants "by the 

values of the modal functions at the stations along the wing and inte- 
grating over the span. The rapidity with which the integrations can he 
performed depends principally on the number of stations needed to 
represent accurately the variation in physical properties along the span 
and on the integration procedure used. 


ACCURACY OF KESUITS 


Because of the structural complexity of airplane wings, simplifying 
assumptions must necessarily he made in any analytical wing -vibration 
analysis. Discrepancies which have been observed between computed modes 
and frequencies and experimentally determined ones will therefore include 
errors due to the assumptions made in the analytical solution and errors 
due to the fact that the analytical solution 1 b usually derived by an 
approximate method. In order that the magnitude of the latter type of 
error might be investigated for the energy solution, a few idealized 
numerical examples are presented for which exact solutions (based an the 
same simplifying assumptions, see appendix C) can also be made. A brief 
description of the examples and the results obtained follows. 

Uncoupled free -free bending vibration . - A simple test of the 
convergence of the energy method for uniform -cantilever-beam mode 
expansions is to calculate the bending modes of a free -free beam with 
the -use of only a few terms to represent the deflection. The following 
deflection expansion 


y = a Q + a^Cx) + a 2 y 2 (x) + a 3 y 3 (x) 


( 15 ) 


leads to a solution for symmetrical free -free modes represented by the 
determinant In table V. The calculation of the constants A-^ , B-j-p, . . . 

for this- case is simplified because a unifoim beam is being used. The 


computations for one of the main diagonal terms 




is, for 


example 
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n l 


2 


my n 


dx 


UO 





dx 


PL 


= m 


jJ 2 - dx 


0 




UO 



p 

where p n Is the nth tending frequency of 
(clamped at the center line) • The integral 


the team as a cantilever 
1L 

p 

y n dx 'has teen 


do 

evaluated (see reference k, p. 336 ), and for tending inodes of a uniform 
cantilever the integral equals L/4 if the modes are given in terms of 
a unit tip deflection- Making use of this fact and dividing through ty 
the factor mL leads to diagonal terms of the form Bhown In tatle V- 
The zeros appear in the tatle tecause of the orthogonality of the modal 
functions in this uniform team example. 


The three lowest frequencies obtained from the fourth -order deter- 
minant are compared with the exact frequencies in the following tatle: 


Symmetrical 

modes 

Frequency 
(radians /sec) 

Percent 

difference 

Exact 

Energy 

1 

5.60 

5.60 

0 

2 

30-25 

30.26 

0 

3 

74.50 

77-46 

4-0 


When the a^y^ term was added to equation (15), a frequency of 
75-37 radians per second ( 1.2 percent error) was obtained for the third 
symmetrical mode- The modes calculated by the energy solution are com- 
pared with the exact modes in figure 1- The solutions for both a 
fourth -order and fifth -order determinant are shown for the third mode- 
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In order "bo calculate the antisymmetrical nodes of the same team, 

Hq in equation (15) is replaced by J-) 'which leadB to the determi- 

nant in table VI. The two lowest frequencies calculated from this 
determinant are compared with the exact solution in the following table! 


Antisymmetrical 

modes 

Frequency 
(radians /sec) 

Percent 

error 

Exact 

Energy 

1 

15.42 

15.42 

0 

2 

50.0 

5O.7 

1.4 


A comparison of the derived modes with the exact modes is shown in 
figure 2. 

From these comparisons the energy solutions are seen to be in good 
agreement with the exact solution for the first and second symmetrical 
modes with a fourth-order determinant and for the third mode with a 
fifth-order determinant (fig. l) • For the first and second antisym- 
metrlcal modes (fig* 2) good agreement is found from a fourth-order 
determinant. In the symmetrical cases, the derived modes contain at 
least one less nodal point than the highest uniform-cantilever -beam mode 
used in the deflection expansion and for the antisymmetrical case, the 
derived modes contain no more nodes than are contained in the uniform 
cantilever modes. (Zero deflection at the root of a cantilever and at 
the center line of an antisymmetrical free -free beam is counted as a 
node.) A simple and, in some cases, conservative rule of thumb, then, 
that might be followed to ensure adequate representation of a desired 
mode containing n nodes is to include in the expansions for deflection 
wl 1 terms up to end including the modal function which has n + 1 nodes . 
This rule should be applied separately to both the deflection and 
rotation components of a mode in the case of coupled vibration. 

Mass coupled vibration of a swept wing and fuselage . - A few of the 
natural frequencies of vibrations of the configuration shown in figure 3 
have been dete rmin ed by both the exact and energy methods. The swept 
wing is of constant croBS section along the span and is mounted on an 
inflexible fuselage. The parameters assumed are as follows: 
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Wing parameters: 

Young *8 modulus of elasticity E, psi 10,000,000 

Shear modulus G, psi 4,000,000 

Bending moment of inertia I, inches^ 800 

Torsion constant J, inches^ 1600 

Mass of wing m, pound-seconds^/inches^ 0.025 

Polar moment of inertia pound -seconds^ l6 


Fuselage parameters: 

One -half of pitching polar moment of inertia iFp, 

pound-seconds^ -inches 

One - h al f of rolling polar moment of inertia Ifr, 

pound-seconds^ -inches 

Distance Between mass center of fuselage and elastic axis 
of wing at root e^, inches 

“jT 


Eatio of fuselage mass to wing mass 


400,000 

10,000 

-70 

3 


An exact analysis of the vibration of this configuration leads to two 
frequency equations, one for symmetrical and one for antisymmetrical 
modes. The derivation of these equations is given in appendix C. 

The energy solution for the symmetrical modes was made "by rise of the 
first four terms in the expansion for deflection (equation (6)) and the 
first three terms in the expansion for twist (equation (7))* Since a Q is 
eliminated in a solution for symmetrical modes, the resulting determi- 
nant, table 711, 1 b of the sixth order. 

Table YU has Been simplified in the manner which has Been previ- 
ously described for tables 7 and 71. In addition, the zeros in the 
upper-right and lower-left quadrants of this determinant are present 
Because there is no mass coupling along the wing. 

The three lowest frequencies satisfying the determinant are compared 
with the exact frequencies in the following table: 
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The modes associated with the frequencies obtained in the energy solu- 
tion are presented in figures k } and 6 * F or the third mode the rule 
of t humb presented previously is not complied with wholly. The per- 
centage error in frequency, however, is rather small hut, in order to 
ensure that the modal deflection is represented accurately, one addi- 
tional term should probably be included in the expansion for deflection. 

In order to check the derivation of the antisymmetrical frequency 
equations the lowest antisymmetrical -mode frequency was found from the 
fourth -order determinant shown in table VTTI. The value was Ik- 66 cycles 
per second which is also the exact value. The associated mode obtained 
from the energy solution is presented in figure 7 * 


DISCUSSION OF METHOD 


As shown by the numerical examples, application of the energy 
method to vibration problems leads to results which compare favorably 
with known exact solutions. Accuracy of results alone, however, would 
not Justify use of the method in a given problem unless the desired 
accuracy of frequencies and modes can be obtained with a reasonable 
amount of labor. Although no attempt is made to compare the labor 
required to compute modes and frequencies by the energy method with that 
required by other procedures, the following facts about the energy 
method and the Crout solutions are presented for consideration: 

(1) The number of terms taken in the expansions for deflection and 
rotation determines the number of coordi na tes in the solution and hence 
the order of the determinants used. The numerical examples have shown 
that solutions for relatively complicated types of wing vibration can 
be obtained from low-order determinants, which is desirable from con- 
sideration of computing time. (Preliminary investigations have shown 
that a number of the coupled modes of the airplane considering the 
fuselage elastic can be obtained from an eighth-order determ i nan t . ) 

Since each mode and frequency is obtained in a manner independent of 
other modes and without recourse to an orthogonality condition, the 
modes may be computed in any desired sequence and with varying order 
determinants , the order depending on the number of terms needed to give 
the particular mode to the deBired accuracy. 

(2) For a given mass and stiffness variation along the wing, most 
of the constants A^, B^, C^, . . . for the wing are used in the 

solutions for both symmetrical and antisymmetrical vibration. Thus a 
given set of equations for symmetrical vibration can be used for 
computing the antisymmetrical modes after only a few new constants have 
been computed. 

( 3 ) The effect of adding to or subtracting from the mas s on the 
wing, of considering elastically mounted motors, or of considering an 
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elastic fuselage is easily introduced into the energy solution. (See 
appendix D for the method of introducing spring mountings. ) 

(^) The Crout method has "been found to he a highly efficient -way 
to solve the characteristic -value equations derived from the energy 
approach. The characteristic values, or natural frequencies, are 
quickly located by a few trial and error evaluations of the determinant 
of the homogeneous equations with assumed values of u in the vicinity 
of a natural frequency; and the unknown coefficients defining the modes 
are then found by a straightforward calculation. 

The expansions for deflection and twist have been made in terms of 
uniform-cantilever -mode functions for several reasons. Other functions 
could have been used. For example, power-series expansions for both 
deflection and twist have been used for the problem considered in this 
paper. These expansions work well, are simple and convenient, but may 
involve computational difficulties which are not present when cantilever 
functions are used. The off-diagonal terns in the determinant for the 
power-series solution are usually larger than the diagonal elements. 
Small differences of large numbers result in the evaluation of the 
deteiminant so that the computations must be carried to a large number 
of significant figures in order to avoid errors. When uniform canti- 
lever functions are used, however, there is in most problems a 
resemblance between each cantilever function and one of the modes being 
derived. The functions therefore have a tendency to be orthogonal 
(since they are orthogonal when the structure is uniform) with the 
consequence that the off-diagonal terms of the determinant are small in 
comparison with the diagonal terms. Evaluation can therefore be made 
easily without the use of a large number of significant figures. The 
examples presented in this paper were computed with the use of four or 
five significant figures. 


CONCEOSIONS 


A solution is presented for the determination of the coupled and 
uncoupled miodes and frequencies of wings (swept or unswept) mounted on 
a fuselage. Characteristic equations describing s ymme trical and anti- 
symmetrical modes of vibration are derived by an energy method which 
makes use of the natural modes of a uniform cantilever beam. The use 
of these modes and the principle that only the geometrical boundary 
conditions need to be satisfied when a solution is made results in a 
direct procedure with a considerable saving in work. The method points 
out that the dynamical boundary conditions, such as shear and moment, 
are automatically satisfied when this principle is used. 
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The application of the characteristic equations to various types of 
•wing Titrations is illustrated "by numerical examples. The numerical 
e x amp les chosen vere susceptible to exact solution and shov that accurate 
modes and frequencies can he obtained from low -order determinants. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., August 30, 19^8 
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APPENDIX A 

f 

DERIVATION OE EXPRESSIONS EOE KINETIC CTERGI 


In the energy solution used herein, the potential energy stored in 
the wing at maximum displacement and the kinetic energy of the wing- 
fuselage system in passing through the equilibrium position must he 
known. Since the equations for potential energy of tending and twist 
are well known (see, for example, reference 4), the present appendix 1 b 
confined to a derivation of the kinetic -energy expression- 

In figure 8, a cross section of the wing is shown at the instant 
that it passes the equilibrium position) the elastic axis is assumed to 
have a vertical velocity v, and the cross section is assumed to be 
rotating at an angular velocity D. Any element of masB dm having the 
coordinates (r,0) can be shown to have a total velocity such that 

v^ 2 = (v + ttr cos 0) 2 + JVrr 2 sln 2 0 (Al) 


The kine tic energy of the element will be 


1 P 

=• dm vt . 

2 


If y and (p are 


the maximum values of deflection and twist, the velocity v and rota- 
tional velcoity 2 may be shown to be equal to coy and o><Pj respec- • 
tively. Substitution of these values in the expression for total ■ 
velocity and integration of the kinetic energy of all the elements over 
the cross section gives for the total kinetic energy of unit length of 
the wing at. the cross section under consideration 


~ m(y 2 + 2eytp + 


where e is the distance between the elastic axis and the center of 
gravity (e Is positive when center of gravity is forward of the 
elastic axis) of the cross section and k is the radius of gyration 
of the cross section about the elastic axis. Integration of the kinetic 
energy over the length of the wing gives for the total kinetic energy of 
the wing 



r 1 

m^y 2 + 2eytp + dx 


(A2) 


u 0 
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Tlie egression for the kinetic energy of vertical and pitching 
motion of the fuselage can "be found "by applying equation (A2) to the 
fuselage ma ss. The angle <p, however, is replaced "by the pitching angle 

of the fuselage given hy ( <p cos A - sin A ) . If 0 is used to 

> /*=0 

denote the pitching angle, the kinetic energy of vertical and pitching 
motion of the fuselage is 

■f + 2e I^ e + k F 202 )] 


or 


p / 

■g'f mpy 2 + SotpepyS + 


' 2 ) 
/x=0 


o 

The kinetic energy of the fuselage in rolling motion is 




where $ is the angle 
^<p sin A + ~ cos A^ 


/ 


then 


x=0 


roll of the fuselage given hy 

The total kinetic energy of tlie fuselage is 


Te-f 


( 


+ 2ntp©jiyQ + Iy 


e 2 + I 




(A3) 


in which y, 0, and $ are understood to he the maximum values of those 
quantities . 
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AEPETTEEX B 


EEEUJITIOn OF COBSTAHTS HJ CHAEACTERISTIC EQUATIONS 


The constants Ain. Binj C in_. . . . found In. the equations of 
tables I and H represent integrals for -which the integrations are 
performed from, the center l in e of the fuselage or root of the wing to 
the ving tip. The integrations are as indicated: 


■Ain = 


®in = 


Cin = 


D in = 


PL 


my^n dr 


0 

PL 


EEyiV n " dz 


UO 

PL 


Im'Pi'Pn 


do 

PL 


GJqj^’qPn’ dr 


do 


B’in ~ 


Gh - 


^n - 


P = 


mejiCPn dr 


do 

PL 


me i dr 

0 

PL 

■m | y n dr 

0 

PL 


m 


(!) e 


dr 


JO 
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APPENDIX C 

HEKEVATION OP FREQUENCY EQUATIONS FOR VIBRATION 
OF SWEPT WINGS 


An exact solution for frequencies and modes of Titration of swept 
wings can readily te made only for an Idea li zed structure representing 
the swepthack-wing airplane. Such a configuration is shown in figure 3 . 
The deflection and twist of the wings are assumed to conform to elemen- 
tary engineering team, theory. The motion of the wings can then te 
derived ty a particular solution to the following familiar differential 
equations 


EE ^ 
dx 4 


(Cl) 


GJ 0 = -I^q) 


(C2) 


The general solutions are 


y = Ci sinh px + C 2 cosh px + C 3 sin px + C 4 cos px (C 3 ) 


tp = Ctj sin qx + Cg cos q_x (C4) 


where 


P 



(C5) 


q. 



(c 6) 


GJ • 
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Values for the sir unknown constants in equations (C3) and (C4) must he 
determined from the known boundary conditions at the root and tip of the 
wing. At the wing tip, the following conditions exist: 


d^y _ d^y d£ _ Q 

dx 2 dx3 dx 


(C7) 


When the wing is vibrating in symmetrical modes, the following conditions 
at the root must also be true: 


/ dy 

V dx 


cos A + <p sin A 


'x=0 


0 


(C8) 



dy 

e — sin A + e<p cos 
dx 



= 0 


(C 9) 



2%stoA + GJ 

dx 2 


^2. cos A + cos A - 


dy 

dx 


sin A 


) 


+ CO' 


fm^epy 


= 0 


x=0 


(CIO) 


Substitution of equations (C 3 ) and (C4) into the six boundary -condition 
equations leads to sir homogeneous equations in terms of the unknown 
coefficients C l> C 2^ * * * and the unknown frequency a>. Since values 

of C-l, C 2 , • • • other than zero which satisfy the set of equations 

can be found only when the determinant of the set of equations is zero, - 
the natural frequencies of vibration are found as those values of co 
which cause the determinant to be zero. Expansion of the determinant of 
the set of equations leads to a single equation containing the unknown 
frequency oj. This frequency equation can be simplified to the form 
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e3i. 


' F P . ™GJ cos 2 A 


+ SB — V - -— tan 0 N 


mL3 eia A n sillA 


(cos 0 sinh 0 + sin 0 cosh. 0) 


4 - (1 + cos 0 cosh. 0)1 - — sin 0 sinh 0 

mL J ml£ 


( p . p 

0 1 + cos g cosll g) + g^Q ^ ( C OS 0 sink 0 

mL 2 / sin A . ' mL 


sin 0 cosh 0) - sin A (1 - cos 0 cosh 0) =0 


(Cll) 


•where p, q, and. <d have "been replaced hy the parameters 0 and B 
having the following relationship to p, q, and cn 


0 = Lp 



(C12) 



(C13) 


Therefore 



5a H 

m QJ 


(C14) 


The roots 0 to equation (Cll) define the frequencies of the natural 
symmetrical modes of vibration of the swept -wing configurations . Bor 
the par ame ters given in the section entitled "Mass coupled vibration- 
‘ of a swept wing and fuselage /' the first three roots to this equation 
are 
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«l = 1*959 
0 2 - 3*336 
©3 = 4.772 


In order to derive the frequency equation for the antisymmetrical 
modes of vibration^ the "boundary conditions defining antisymmetrical 
Titration must "be satisfied. In addition to the "boundary conditions 
given "by equation (C7), the following conditions must "be true: 


(y)x=o = 0 


(C15) 


0 


cos A 


_ 

dx 


sin A 


)x=0 


= 0 


(Cl6) 



cos A + GT ££ sin A 
dx 


+ ^ ~ cos A + 9 sin A 


) 


= 0 


x=0 


(CI7) 


Use of the conditions given in equations (C7), (CI5), (Cl6), and 
(CI7) yields as the frequency equation for antisymmetrical modes 


— 2 — + 0IT — tan 0%Vl + cos 0 cosh 0) 

,mL3 COB A EE cos A 


+ cos A (sin 9 cosh 9 - cos 9 slnh 0) = 0 


(Cl8) 
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For the parameters given in the section entitled "Mass coupled vibration 
of a swept wing and fuselage," the first three roots to this equation are 


©1 = 2.551 


e 2 = 3.462 


e 3 = 4.79 


0 
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appendix d 

THE EFFECT 03? A SPRING SUSPENSION ON THE NATURAL FREQUENCIES 
AND MDEES OF A WXNG (SWEPT OR UNSWEPT) AND FUSELAGE 


Frequencies nnfl modes of airplane •wings are often deteimined experi- 
mentally with the airplane suspended at the center of gravity hy a linear 
spring. An evaluation of the changes in vibration characteristics 
(symmetrical modes only are affected) of the airplane due to the spring 
suspension can easily he made hy the energy method. The potential energy 
stored in the spring is given hy 


a ( y + 


e a ( P 

cos Ay 


x=0 


(EL) 


With the expansion used for y and <p, equation (Dl) reduces simply to 




e gkp \ 
cos A / 


(D2) 


If Y a is added to the energies U, Tq_, and Yq } and the expression 
U + Y a - Y]_ - Yq is minimized with respect to the coefficients a Q 
and h Q , the following equations are derived: 

For a Q , 


®o/W 








mj©ji 

cos A 


- LHj_ tan A 


V t 

- — ) + Y” hnFia = 0 (i = o) (»3) 

cd2 cob A J 
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For 


J o’ 


' g 

a o /F 10 + - LHi tan A - — -5SelJ\ + ]T a^F^ - 1^ tan a) 

\ cos A <xr cos A/ n=l ' ' 


/ 2_2 ^Fp 06, 

+ T) 0 I w + LT? tan^A + r 


cos^A 


2TGi tan A - 


_j£J) 

03^ cos^Ay 


+ ' b n( c 3n - 1% tan A) = 0 (i = 0) 


(B4) 


These two equations replace rows for ao and b 0 in table I when the 
natural vibration is restrained by a linear spring at the center of 
gravity of the airplane. 



28 


NACA TFT No. 17^7 


1. Lawrence, H. R.: The Dynamics of a Swept Wing. Jour. Aero. Sci., 

vol. 14, no. 11, Nov. 1947, pp. 64-3-650. 

2. Courant, R. : Yariatlonal Methods for the Solution of Problems of 

Equilibrium and Vibrations. Bull. Am. Math. Soc., yoI. 4-9, 

Jan. 194-3, PP* 1-23* 

3. Crout, Prescott D. : A Short Method for Evaluating Determinants and 

Solving Systems of Linear Equations with Real or Complex Coef- 
ficients. Trans. A.I.E.E., vol. 60, 194-1, pp* 1235-1240. 
(Abridged as Merchant Methods MM-182, Sept. 194-1, Merchant 
Calculating Machine Co., Oakland, Calif.) 

4-. Timoshenko, S.: Vibration Problems in Engineering. Second ed., 

D. Van Nostrand Co., Inc., 1937* 



xjtt&z i.- jqotnajs it* ssi«B 2 hidal ttbeahoh » met ran 


jror n± | 

i - ° faw + ”?) ^oi *oa •'cc ^bo + ” 130 *“ A ) *01 *oa ■ r ° m 

l r. l j >oi \Au. - (ale - -^) (Ain * ~^}J ( T io - A ) J i» 

l 

f 

1 

i ■ a | *oe (*12 ‘ (hz ' (a 2n * "jr) C'eo - ^ t<m $ 'ai T s2 ? e* 

( _> 

1 “ n Ajn ^ln * -^) (Ate - (\c ' -^) (*00 - ten a) T^i r n2 'mi 

I 

»or bi 

i U 0 ^?bo + - IH 0 ten A^ (rio -I^tuj) (Ted - Hfe tm A) ^nO - £% tan A) ^0<x> + L 0 * ten?A + - filbo tan A^ (0 01 - 1!^ tn i) C°oe * LOe tan a) (0^ - Ifl>, tan a) 

1.1 r 01 in Ita. (°oi - ifli tan a) (o u - ^ (<>i£ - ^-f) (<>i* - ^f) 

i » e Ha Tib the Hi2 (Poa ■ iba ten a) (oje - (^sa • (® 2 « ■ 

1 . ■ Tqc Tl* Tea Tom (°Cn - Iflc tan 4 ) ^o» * “^) (°2e ‘ ^f) (“nm ■ 


t 


LU I 'Oil iE> VOVK 



TABLE H. - KSMTIOIB FOR AKTISTWffiTRICAL VXBHJCTXOB 0? 3UEPT «U» 


LO 

O 



LijLT 'OK Mj VO VH 



TABLE III. - BENDING- MOIES OP A UNIFORM CANTILEVER BEAM AND THEIR SECOND DERIVATIVES 


X 

L 

n 

72 

y 3 

n 

y 5 


1.0 

1. 0000 

1.0000 

1.0000 

1.0000 

1.0000 

0 

•9 

.8624 

.5237 

.2285 

-.0520 

-.2940 

1 

.8 

•7255 

.0700 

— 39^9 

-.6430 

-.6005 

.2 

*7 

.5909 

-3170 

“•65 7 ^ 

- 397 ^ 

.2257 

•3 

.6 

.4611 

-.5895 

-.4738 

.3265 

.7002 

.4 

•5 

•3395 

-7136 

.0197 

.7071 

.0007 

•5 

.4 

.2299 

-.6835 

•5297 

•3156 

-.6966 

.6 

• 3 

.1365 

-.5261 

•7562 

-4339 

-2113 

hqh 

.2 

.0638 

-.3030 

.6045 

-7538 

.6596 

.8 

.1 

.0168 

-.0926 

.2281 

-3850 

.5372 

•9 . 

0 

0 

0 

0 

0 

0 

1.0 


I 2 *,; 

3*516 

lV 

22.03 

L %'' ■ 
61.70 

_ iV 

120.9 

lftr 5 " 

199.9 

1 

L 
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TABLE V.- 


*0 


1 

1 L 

7i cLz 

Jo 

TL 

y 2 dx 

Jo 

">L 

y 3 ax 

Jo 


p x 2 = 12 
TL 

1 y x dx = 
UO 


EQUATIONS FOE SYMMETEICAL FEEE -FREE -BEAM EXAMPLE 


1 

L 


1 

4 


a l 

®2 

1 L 

1 L 


y 2 dx i 

0 1 

0 l 

i-*£) 

0 




PL 


y 3 


dx 


0 



0 


0 


0 





36 P 2 2 = 485.5 P3 2 = 3807.O 


•3915 


PL 



y 2 ax 


-O.2169 




tr'IH 


3 ^ 
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TABLE 71.- EQUATIONS FOE ANTISIMMEIRICAL EKEE -FREE -BEAM EXAMPLE 


a l «2 a 3 



i 


' o 
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TEffTP. 7H. - EQUATIOIB FOR flTK gg P RTO A T, YTBRATim OF SHZPT-WDK3- EIAMPIB 


®o 


a l 

»e 




1L 

PL 


L + iy 

1 

L 

1 

Jl dx 

0 

*1*“ 

t 

'L 

Jl dx 

0 

K i -S) 

0 

. 1 

’L 

7 2 41 

0 


0 

}(*■« 

- 


\coe A 2/ 


PI 


L^l' 


PL 




(V£_l\ 


\ooo A 2 / 

l 


PL 


S'i fc 


PL 


L^ 41 


fi + i£ + _^_ N ) H 

\ “■ 3 m y ooo 2 A j * L 


*i< 


ii 

a L 


PL 




di 


do 


il 

m L 


dx 


8 » V a?/ 


Jsi. 

m L 


92 


ifefi .2^ 

2 » V 


p, 2 - 12.36 - e472.h- 

1 mL 4 


^i 2 " f “2 - 2 ’ 00 ” 2 

tt 


ft, 2 . bB6 K - 97105 

mL 4 




^'S“ 3 


Js . 6to 






Ao 45 0 


tan A- 1 


(Bo© table* Y and YI for Tallies of other Integrals.) 



TABLE M. - EQUATIONS IDE ANTISI1METBICAL 7IBEATI0N OF A SWEPT -WING EXAMPLE 


U) 

o\ 


% 


a l 


“a 


+i 


i (x'f ax + i + Fr 

L \L/ Kw-L cos^A, 

OO 


1 

L 


E r l 41 


JO 


n t 


1 

L 


|y 2 to 


ii 

m L 


n t 


<Pj_ dx 


1 

L 




l(l „!£ 

o£ 


1 

L 


' Jj 


— yw dx 

L 




Vi f 


m L 


*1 


dx 


0 


1 . %L C 

2 a \ a£ 


P l 2 = 2k2TJ.k 


p 2 2 = 97105 


% 


2 B IQOOOit 2 

k 


“o' 


tea A 


A = 45 c 


(See tables YI and YU for values of Integrals.) 
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Figure 1.- Symmetrical modes of uniform free-free beam 


Relative deflection Relative deflection 
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Figure 2.- Anti symmetrical modes of uniform free- free beam. 
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Figure 5»- Example swept-wing airplane 
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Figure First symmetrical mode of example swept wing. 



Figure 5 .- Second symmetrical mode of example swept wing. 
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Figure 6 .- Third syranetrical mode of example swept wing. 



Figure 7 .- First anti syranetrical mode of example swept wing. 
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Elastic axis 



Figure 8.- Coordinate and velocity notation for 
element of mass on cross section. 



